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Exciton fine structures in cubic III-V semiconductor GaAs, InAs and GaN quantum dots are 
investigated systematically and the exciton spin relaxation in GaN quantum dots is calculated by first 
setting up the effective exciton Hamiltonian. The electron-hole exchange interaction Hamiltonian, 
which consists of the long- and short-range parts, is derived within the effective- mass approximation 
by taking into account the conduction, heavy- and light-hole bands, and especially the split-off 
band. The scheme applied in this paper allows the description of excitons in both the strong and 
weak confinement regimes. The importance of treating the direct electron-hole Coulomb interaction 
unperturbatively is demonstrated. We show in our calculation that the light-hole and split-off bands 
are negligible when considering the exciton fine structure, even for GaN quantum dots, and the short- 
range exchange interaction is irrelevant when considering the optically active doublet splitting. We 
point out that the long-range exchange interaction, which is neglected in many previous works, 
contributes to the energy splitting between the bright and dark states, together with the short- 
range exchange interaction. Strong dependence of the optically active doublet splitting on the 
anisotropy of dot shape is reported. Large doublet splittings up to 600 j-ieV, and even up to several 
meV for small dot size with large anisotropy, are shown in GaN quantum dots. The spin relaxation 
between the lowest two optically active exciton states in GaN quantum dots is calculated, showing 
a strong dependence on the dot anisotropy. Long exciton spin relaxation time is reported in GaN 
quantum dots. These findings are in good agreement with the experimental results. 
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I. INTRODUCTION 

Semiconductor quantum dots (QDs) have attracted in- 
tense interest due to the high potential to work as basic 
device units for photonics, spintronics, quantum commu- 
nication and computation^— One of the applications is 
to use QDs as emitter of single photon 8 - - — and entan- 
gled photon pairs 1 ^— via the combination of exciton and 
biexciton. Explicitly, for a symmetric QD, the decay of 
the bright (optically active) exciton states^ 9 - with total 
angular momentum projections J z = ±1 can emit a 
circularly polarized photons, while that of the biexciton 
states through two intermediate degenerate bright exci- 
ton states can generate two entangled photons i 12 ' 20 The 
asymmetry of QDs lifts the degeneracy of the two bright 
exciton states and mixes them into two states that gener- 
ate linearly polarized photons with orthogonal polariza- 
tions during the decayi 15 ' 16 : 21 " — This is even crucial when 
considering the biexciton decay as the splitting of the in- 
termediate exciton states makes the two channels distin- 
guishable and hence destroys the ent anglement 1 14 i 18 i 23 i 24 
Moreover, the dark (optically forbidden) exciton states 
which lie slightly below the bright ones 1 ^ are consid- 
ered promising candidates for spin qubits due to their 
long lifetimes when confined in QDs i 25 ' 26 They also 
play a key role for the Bose-Einstein condensation in 
semiconductors £1 

The splitting between the bright and dark exci- 
ton states (denoted as the BD exchange splitting 



hereafter) 28 ' 29 and that between the bright exciton 
states (denoted as the doublet splitting hereafter )22r— 
are mainly controlled by the electron-hole (e-h) ex- 
change interaction together with the dot size and 
asymmetry j 22 i 32 i 36 ~ — So a detailed understanding of the 
e-h exchange interaction and the resulting exciton fine 
structure in QDs is of fundamental importance for both 
theoretical and application purposes. 

The e-h exchange interaction has been investigated 
ever since the 1960s. It can be decomposed into long- 
and short-range parts in the real space 4 ^— or the ana- 
lytical and nonanalytical parts in the k space d 3 ' 44 There 
is a close correspondence between them which can be 
found in Ref. l45l and sometimes no difference is made 
between these two approaches 4& 

Early investigations on the e-h exchange interaction 
mainly focus on the bulk system] 43 ' 47 " — It is well known 
that, when considering the exciton states in semiconduc- 
tors by taking into account the conduction band r§ and 
the valence band Tg, the short-range (SR) exchange inter- 
action splits the eightfold-degenerate exciton state into a 
triplet bright state and a quintuplet dark state with the 
splitting energy between them, the so-called exchange 
energ y 52 ' 53 The long-range (LR) exchange interaction 
further splits the triplet into a longitudinal and two trans- 
verse modes, the energy difference of which is denoted 
as the longitudinal-transverse splitting^ 2 .^ The e-h ex- 
change interaction, as well as the direct Coulomb inter- 
action, is greatly enhanced by the quantum confinements 
in low-dimensional semiconductor structures because of 
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the increased spatial overlap between the electron and 
hole wave functions. Reexamination of the e-h exchange 
interaction in low-dimensional structures was intrigued 
in both experimental^— and theoretica l 19 ' 37 ' 60 " — ways. 
Most of the theoretical works were carried out within 
the framework of the envelope-function approximation 
together with the effective-mass approximation] 37 ' 61- — 
In general, due to the different effective masses of the 
heavy, light and split-off holes, the heavy-, light- and 
split-off-hole exciton states 6 ^ are energetically split when 
a quantum confinement is applied. For common cubic 
III-V semiconductors, the heavy-hole exciton which en- 
ergetically lies the lowest is of most physical interest and 
is hence mostly investigated^ 6 -^- 6 ^^^- 7 ^ The heavy- 
hole exciton quartet, which is characterized by the total 
angular momentum projections J z = ±1, ±2, is split 
into bright and dark exciton states with J z = ±1 and 
J z = ±2, respectively, by the e-h exchange interaction 
and the bright exciton states | ± 1) are further split un- 
der anisotropic confinement potential. Chen et al£^ cal- 
culated the exchange energy and the BD exchange split- 
ting in GaAs/Gai-zAl-rAs quantum wells (QWs) and ev- 
idenced the enhancement of the exchange effect with de- 
creasing well width. Their work was based on the approx- 
imation of decoupled heavy- and light-hole subbands. 
Andrcani and Bassani investigated the k dependence of 
the e-h exchange interaction in QW systems^ In Ref.l63l 
the exciton spin dynamics in GaAs QWs was studied with 
e-h exchange interaction as an effective spin-flip mech- 
anism, where the matrix elements of the LR exchange 
interaction were calculated by using the simple heavy- 
hole exciton ground states while those of the SR ex- 
change interaction were obtained by taking into account 
the mixing of heavy- and light-hole bands. Takagahara 
performed systematic studies of exciton states in QDs in 
Ref. [H. The exciton binding energy and the LR and 
SR exchange interactions were investigated using wave 
functions calculated by the variational method. The sub- 
band mixing induced by the direct Coulomb interaction 
was pointed out to be important. In a later work for 
GaAs QDs* 3 - 7 - Takagahara derived an eight-band exciton 
Hamiltonian. The LR exchange interaction, which was 
attributed to dipole-dipole interaction, was emphasized 
to be much more important than the SR exchange inter- 
action when concerning the exciton doublet fine struc- 
ture. The contradiction of the calculated scaling law of 
the doublet splitting energy to the theoretical prediction 
was reported (also in Ref. [69t ). Whereas in the latter 
works by Tsitsishvili et a?j 66 ' 72 the exciton spin relax- 
ation in single asymmetrical QD was studied by taking 
into account only the SR exchange interaction. Efros et 
al£^ and Horodyska et al. in their very recent work- 6 - 7 - in- 
vestigated the band-edge exciton states in spherical QDs 
by including only the SR exchange interaction. But actu- 
ally, as will be shown in this paper, the LR exchange in- 
teraction contributes to the splitting between bright and 
dark exciton states even in isotropic QDs where the dou- 
blet splitting energy vanishes. Its contribution to the BD 



exchange splitting can be comparable with that from the 
SR exchange interaction. Moreover, although the direct 
Coulomb interaction is believed to lead to poor conver- 
gence if treated perturbatively when the exciton system 
is in the weak confinement regime, in a very recent work 
by Kadantsev and Hawrylak^ the exciton fine struc- 
ture in GaAs QDs was still studied by taking the direct 
Coulomb interaction as perturbation. 

Thus, even though a lot of works have been done on 
the e-h exchange interaction induced exciton properties, 
obvious confusions are still seen in the literature, even for 
the most investigated Ini-^Ga^As nanostructures. The 
significance of valence-band coupling to the exciton fine 
structure needs to be evaluated. The relative importance 
of the LR and SR exchange interactions to the doublet 
splitting needs more elaboration, and that to the BD 
exchange splitting needs to be clarified. Furthermore, 
the size scaling of the e-h exchange interaction and the 
influence of the direct Coulomb interaction on the exciton 
fine structure have to be examined and stressed. 

Works discussed above mainly concern III-V semicon- 
ductor nanostructures based on Ini-^Ga^As the split-off 
band of which is far away from the heavy- and light-hole 
bands ^ So the split-off band is always neglected when 
studying exciton fine structure. However, for cubic GaN, 
whose spin-orbit splitting is small compared to the wide 
band gap? 7 - 3 - no theoretical work has been performed to 
investigate exciton fine structure so far, nor does the ex- 
plicit expression of e-h exchange Hamiltonian with the 
effect of split-off band exists in the literature. It has been 
proved that in bulk GaN, the split-off band is important 
when considering the spin-orbit coupling. 74 Whether it is 
still the case when studying the e-h exchange interaction 
in GaN QDs needs to be examined. 

The exciton spin relaxation is another important sub- 
ject that strongly affects the quality of applications in in- 
formation storage and processing based on exciton states 
in QDs i 66 ' 75 " — Motivated by recent experimental study 
of exciton spin orientation in cubic GaN/AIN QDsjS we 
investigate the behavior of the LR and SR exchange inter- 
actions in cubic III-V semiconductor QDs, taking into ac- 
count all heavy-hole, light-hole and split-off bands. The 
spin relaxation between the lowest two bright exciton 
states in single GaN QD is studied after that. 

In this paper, in order to explicitly include the effect 
of the split-off band on the exciton fine structure in cu- 
bic GaN, GaAs and In As QDs, the bulk e-h exchange 
Hamiltonian of both LR and SR parts is first derived 
in the 12 x 12 matrix representation. The derivation is 
carried out within the framework of the effective-mass 
approximation. 79 For systems strongly confined in one 
direction, e.g., the QD system with small dot height con- 
sidered in this paper, we are able to apply the Lowdin 
partitioning method 7 - 9 -"— to approximately diagonalizc 
the modified 6x6 Luttinger Hamiltonian (specified in 
Sec. II) for holes to obtain a new "heavy-hole" band, 
which is an admixture of the heavy-hole, light-hole and 
split-off bands. In this way, 4x4 matrix representations 
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of the exciton exchange Hamiltonian are constructed by 
taking the conduction band and the new "heavy-hole" 
band with the effect of valence-band mixing included. 
We then apply the effective Hamiltonian obtained to in- 
vestigate the exciton fine structures in cubic III-V semi- 
conductor QDs. The doublet splitting energy and the BD 
exchange splitting are calculated and the relative impor- 
tance of the LR and SR exchange interactions as well as 
that of the heavy-hole, light-hole and split-off bands are 
discussed. The size scalings of the doublet splitting en- 
ergy and the BD exchange splitting contributed by the 
LR and the SR exchange interactions are analyzed and 
explained by the scaling rules established. 

Due to the fact that the values of the exciton Bohr ra- 
dius in bulk GaAs, InAs and GaN are f4.9 nm, 51.6 nm 
and 4.8 nm,— respectively, which are comparable or 
even smaller than the average diameter of the QDs in 
this paper (given in Sec. Ill), chosen according to the 
experiments , 29 ' 38 ' 78 ' 83 the direct Coulomb interaction is 
too large to be treated perturbatively. We solve the 
Schrodinger equation by taking the direct Coulomb in- 
teraction and the confinement in an equal footing. The 
calculated exciton binding energies are found to be con- 
siderably large, markedly enhanced by the confinement. 
This contradicts the results in the latest work by Kadant- 
sev and Hawrylak 69 and further demonstrates the impor- 
tance to treat the direct Coulomb interaction unpertur- 
batively. The importance of the direct Coulomb inter- 
action to the exciton fine structure is demonstrated. Fi- 
nally, the exciton spin relaxation rates in single GaN QDs 
are calculated and long spin relaxation time is obtained. 
Our results are in agreement with experiments i 31 ' 38 ' 78 

This paper is organized as follows: In Sec. II, we set up 
our model and lay out the formalism. Matrix represen- 
tations of LR and SR exchange interactions are derived 
first in bulk and then in QW. Size-scaling rules are estab- 
lished and exciton spin relaxation assisted by the acoustic 
phonons are introduced after that. The numerical scheme 
is laid out at the end of this section. In Sec. Ill, we show 
our numerical results of exciton fine structures in GaAs, 
InAs and GaN QDs. Explicit properties of the LR and SR 
exchange interactions contributing to the doublet split- 
ting energy and BD exchange splitting are discussed in 
detail. In Sec. IV, the exciton spin relaxation in single 
GaN QD is investigated. We conclude in Sec. V. 



cal scheme is laid out at the end of this section. 



A. e-h exchange interaction 

1. e-h exchange interaction in bulk 

We start our investigation on the e-h exchange inter- 
action from a general description of direct Wannier-Mott 
excitons in bulk system within the framework of effective- 
mass approximation. The exciton wave function can be 
written in the for m 41 i 42 

*(ri,r 2 ) = ^ [fmri(ri,r 2 )Vw (ri)?/'rik (r2) 

inn 

(r 2 ,ri)^ roko (r 2 )V! n ko(ri)], (1) 

where V'mk (r) is the conduction-band Bloch function 
and ip n k ( r ) 18 the Bloch function for the hole which 
is the time reversal of the Bloch function of the miss- 
ing electron ! 41 ' 42 As for cubic III-V semiconductors, we 
are interested in excitons at the band edge, i.e., the T 
point with ko = 0. m (n) is the index for the electron 
(hole) band under consideration, including the spin de- 
gree of freedom. F mn {Yi,Y2) is the envelope function 
with -F m n (i"i, 1*2) = —F nm (r 2 , i"i). This makes the exci- 
ton wave function antisymmetric. 

The eigen equation for the envelope function F mn is 
given b y 41 ' 42 



E 



dridr 2 ffl h , 
= EF m , n ,(r[,r' 2 ) 



r l r 2 

ri r 2 



F mn {r 1 ,r 2 ) 



(2) 



The explicit form of H c ^, n , [Z\ T T \\ is given in Ap- 
pendix A. 

Now we proceed to a more detailed derivation of the 
matrix representations of the exchange interaction for 
cubic III-V semiconductors, such as GaAs, InAs and 
GaN, by taking into account the conduction band Tg, 
the heavy-hole and light-hole bands Fg and the split-off 
band Tj. The Bloch functions for these bands and the 
time reversal of those for the valence bands are given in 
Appendix BM- We further denote the conduction band 
Bloch functions as 



II. MODEL AND FORMALISM 

We formulate the theory of excitons in cubic III-V 
semiconductors by taking into account the conduction 
band T§ and the valence bands Tg and Ty. The effec- 
tive representations of the exciton Hamiltonian for the 
LR and SR exchange interactions are derived first in 
bulk and then in QW. The size-scaling rules of the e- 
h exchange interaction are established after that. We 
then introduce the exciton spin relaxation due to the 
electron/hole-acoustic-phonon scattering. The numeri- 
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'2' 2> 



|c 2 ) = l|,-|> 



and the hole Bloch functions as: 
, ,3 3, , , ,3 1, 



\ v i) = \^ + ^)h, lV2) 



\2> + 2 )h, lV3) 



(3) 



3 3 11 11 

K) = \^,-^)h, \v 5 ) = \^,+^)h, k) = \^,-^)h- (5) 

We perform the Fourier expansion as U(r) = 
J-x f dqC/ e lq r where U a = e ■> , and then the term 
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of LR exchange Hamiltonian in Eq. (|A6[) transforms into 
the form 

a/3 

x e^(n-ri) ( 5( ri _ r^gft „ ^ ( 6 ) 

with being the time reversal operator. We define the 
band- relevant part of the LR exchange Hamiltonian as 

Qm'n' (q) With 

Q m ^(q)=^Q^e»fe^- (7) 

a/3 e "' m 

I 



Then the parity of the Bloch functions enables us to 
write the matrix representation of Q m i„i (q) in a some- 
what simple way. With the basis \mn) taken in the order 
\cxVx), \c 2 vi), |ciw 2 ), \c 2 v 2 ), \ciV 3 ), \c 2 v 3 ), \civ 4 ), \C2V4,), 
\civ 5 ), \c2V5}, \civ e ), \c 2 v 6 ), Qm'n' (q) takes the form 
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The other half of the matrix is obtained by taking the 
Hermitian conjugate. The matrices below are given in 
the same way. Here K± = q x ± iq y and K = {q x ,%)- 
For the coefficients, A = n h fL , B = , 2 = 'fj' 2 , . — r 

and C = Tnk ^ A — ^2 with too, E„ and A so stand- 
ing for the free electron mass, the band gap and the 
spin-orbit splitting, respectively^ P = (S\p x \X) = 



(S \p y \ Y) = (S\p z \Z). ?£1= E P with E P being a band 
structure parameter in energy unitj^ Approximation has 
been made that the element R = 4r ^ e2 (S \-^Vq \ X) = 

orbit coupling in ir [given in Eq 
it is always much smaller than P. The 12 x 12 matrix 
representation of SR exchange interaction is written as 



(S\£V \Z) from the spin- 
4}], is neglected since 
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in which 

D = V^ I ^K^ q ' r P0l 2 . (io) 

An 8x8 matrix representation of the LR and SR ex- 
change interactions can be found in Refs. |63| and l37t 
where only the heavy- and light-hole bands are included. 
They are the same as the 8x8 submatrices at the top 
left corner of our expressions, but in two— or zero- 
dimensional^ forms. Whereas in Ref. [f^l a 4 x 4 exciton 
Hamiltonian was presented including only the heavy-hole 
band. So the previous expressions correspond to only 
part of our results, where the effects from the split-off 
band and even the light-hole band are absent. In compar- 
ison with the matrix form of exchange interaction Hamil- 
tonian in Ref. [6^, we obtain 

2e 2 h 2 E P , s 

37re Km a^ ohr £;^ 

D = J™| ohr A£ SR , (12) 

where E g is the band gap, ASlt and A£sr are the 
longitudinal-transverse and the singlet-triplet splittings 
in bulk which can be obtained by experiment . 62 1 82 m* 
is the effective mass of the conduction electron, 71 is 
the band parameter and aBohr = 47reoK?i 2 /i/e 2 with 
/.t -1 = l/m*+7i/mo is the exciton Bohr radius in bulk. 82 

2. e-h exchange interaction in QW 

With the e-h exchange Hamiltonian described above, 
one can investigate the properties of excitons by explicitly 
including the contributions of the heavy-hole, light-hole 
and split-off bands. In QWs with small well width, as 
an usual procedure, two-dimensional (2D) exciton bound 



states are used to approximately count the effect of direct 
Coulomb interaction and the e-h exchange interaction 
is treated perturbatively . 61 ' 63 For the QD system in the 
strong-confinement regime, in the literature, people first 
solve the confinement and then treat the direct Coulomb 
and e-h exchange interactions as perturbation.-?' 70 How- 
ever, in QWs or QDs the characteristic size of which 
is comparable or even larger than the exciton Bohr ra- 
dius ciBohr, which means that the Coulomb interaction 
tends to overtake the confinement, one has to solve the 
Schrodinger equation with both the confinement poten- 
tial and the direct Coulomb interaction included. In 
this case, it can be extremely CPU expensive to employ 
the 12x12 e-h exchange Hamiltonian. So the Lowdin 
partitioning methoct^r— is employed to derive the e-h 
exchange interaction in a smaller Hilbert space, while 
still taking into account the confinement-induced valence 
band mixing. 

We start from a system with strong confinement along 
the z direction (i.e., the [001] direction). The infinite 
square well potential is employed. The QW width is de- 
noted as l z . For the envelope function in the z direction, 
only the lowest subband is relevant for both electron and 
hole. 

It is noted that the hole Hamiltonian is not described 
by the so-called Luttinger Hamiltonian itself^ which is 
the Hamiltonian for valence electrons^ One has to trans- 
form it to hole space according to rules given in Refs. |4l| 
and The obtained hole Hamiltonian in the valence 
bands Ty and Tg takes the form as —1 times the 6x6 
Luttinger Hamiltonian . 79 ' 85 After applying a strong con- 
finement along the z direction, the Luttinger Hamilto- 
nian is deduced to 2D form where the odd terms of k z 
vanish and (k 2 ) = tt 2 /I 2 . The values of the minima of 
the diagonal elements for the heavy, light and split-off 
holes are now separated due to their different effective 
masses in the z direction^ The heavy-hole band lies en- 
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ergetically much lower than the other two. This enables 
us to apply the Lowdin partitioning 7 - - — and get decou- 
pled new basis functions for holes. The lowest subbands 
are heavy-hole-like states, which are admixtures of the 
heavy- with the light- and split-off-hole states. 

The Lowdin transformation of the hole Hamiltonian is 
given by H h = e~ s H h e s where S is an anti-Hermitian 
6x6 matrix. The basis functions transform as ip n = 
^] m (e )mnipm- Up to the first-order approximation, one 
obtains ip n = X) m (<W + S mn )tp m , with 

&mn being the 

Kronecker delta. Up to the second order, the effective 
Hamiltonian of the new heavy-hole-like subbands takes 
the form 



H 



HH 



h 

h 



(13) 



h 2 h 2 TT 2 

{Ji+J2)k 2 + - 55(71-272), (14) 



2m 



with fcj| = (k x ,k y ), k x — —id x and k y = —id y . The 
non-zero elements of the matrix S up to the first order 
read 



for the energy splitting between the heavy-hole and light- 
hole (split-off) subbands. The other half of S is obtained 
from the relation — — S. New Bloch functions for the 
heavy-hole-like states become 



\Vl) = \vi) + S 31 \v 3 ) + S 61 \v 6 ) 

\V 2 ) = \v 4 ) + S 2i \v 2 ) + S M \v 5 ) 



(21) 
(22) 



Since the major component of \Vi) {\V 2 )) is (|«4)), we 
still denote the coupled states \Vi) and \V 2 ) as spin ±| 
states for simplicity in the following. The new "heavy- 
hole" bands are now admixtures of the heavy-hole, light- 
hole and split-off bands. It is this band-mixing effect that 
makes the dark exciton states, constructed as |ciVi) or 
IC2V2), become partially optically allowed z&gL&l&s 

Now we are ready to derive the e-h exchange interac- 
tion in 4x4 matrix representation with |ctV^*) as the new 
Bloch wave function from Eqs. © and ([5]). In order to 
write the Hamiltonian in a simple way, we note that since 

(x'y',n zl \A\xy,n z2 ) = {n zl \z')(x' y' , z'\A\xy, z)(z\n z2 ) , 
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where A is an arbitrary operator and the Einstein sum- 
mation convention is presumed, we can always write the 
exchange Hamiltonian in three-dimensional (3D) space 
while keeping in mind that the formulae hold true only 
for QW or QD systems with strong confinement in one 
direction. 

The LR exchange interaction, given in the basis taken 
in the order |ciVi), IC2F1), IC1V2), IC2V2) (or expressed as 
the eigenstates of J z : | + 2), | + 1), | — 1), | —2)), is written 
as 
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Here 7 j are the band parameters^ and E ao (E ac ) stands in which 
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"61 V 


V V6 




-r' 2 ), 






SeiY 






r 2 ), 








V V6 





(24) 



(25) 



(26) 



Here £*■ is the complex conjugate of Sy denned in 
Eqs. (|T5| - (|18p with k — > k 2 , and Sy is the same as Sij 
but with k — > k 2 . The matrix of short-range exchange 
interaction is given in the same way: 



h3„ 



in which 



ri r 2 





H 22 







L>f 3 * 



(27) 



0/ 



SR 



iL 



o-SR 
H 23 



-V2D 



H, 



SB 
33 



^(n-r^n-r'O^ra-r^), 

o/ o* 
1 _ J 61 , °24 , °54 

VH V3 \/6 \/3 
x 5(n - r 2 )J(ri - ri)<5(r 2 - r 2 ) 

D*(ri - r a )<y(n - r[)6{r 2 - r' 2 ). 



(28) 

(29) 
(30) 



From the above expressions, one can see that the four- 
fold-degenerate exciton states are split when the e-h ex- 
change interaction is taken into account. When the diag- 
onal matrix elements are included, the quadruplet is split 
into two doublets: the bright and dark doublets with the 
bright one lying above the dark one. The off-diagonal 
matrix elements further couple the two bright exciton 
states together and cause the doublet splitting. More- 
over, it is noted that in Eqs. ([23 ]) -([2l) ]l and ([2"8 > ([3"0 j) the 
terms without S*j and S[j are derived from the heavy- 
hole band; the terms with S^ 2 4) and S' 31 ^ 2 4) are derived 
from the light- hole band and the terms with Sg^^ and 
^61(54) are derived from the split-off band. It is noted 
that, as shown in Eq. (|27p. the SR exchange interaction 
now directly couples the J z = ±1 exciton states due to 
the confinement-induced band mixing. This coupling is 
missing when only the heavy-hole band is included as in 
Ref. M 

We note that the exciton wave function ^(ri,!^) in 
Eq. (fl]) is actually the same as Eq. (2.1) in Ref. [33 up 



to a conventional constant, while the latter is written in 
the framework of second-quantization. This brings some 
insights into how the electron and hole share the prop- 
erties of identical particles. The exciton wave function 
Eq. (fTJ) is also widely used as a truncated on o 61 ' 65 ' 82 ' 89 

*(r e ,r h ) = F cv (r e , r h )tp cko (r e )^ ko (r h ), (31) 



which is sufficient to describe most of the properties of 
exciton, e.g., spin dynamics 6 ^ and fine structure! 61 ' 65 We 
use this expression in the calculation of the exciton spin 
relaxation rate. 



B. Scaling of the e-h exchange interaction 

The study of size scaling of the e-h exchange interac- 
tion in semiconductor QDs helps to understand the ex- 
perimental result s 28 ' 57 ' 90 and provides an intuitive under- 
standing of how it varies with dot sizei 37 ' 58 ' 65 ' 69 ' 91 " — The 
scaling rules were established in the strong confinement 
regime ; 37 ' 65,69 ' 89 For the doublet splitting, Takagahara^ 7 . 
and Kadantsev and Hawrylak 6 ^ established that the 
LR exchange interaction which determines the splitting 
scales as 1/L 3 with L standing for the characteristic size 
of the QD. Their numerical results (fit to C/L n with 
n«1.3 in Ref. [13 and n = 1.3~1.5 in Ref. [691 ) showed 
discrepancy from the 1/L 3 dependence and the discrep- 
ancy was attributed to the details of the envelope func- 
tions. For the BD exchange splitting, most works were 
carried out retaining only the SR exchange interaction 
and therefore the BD splitting was assumed to scale as 
1/L 3 i 28 ' 65 ' 89 ' 91 However, clear deviation of the experi- 
mental results from the 1/L 3 law was reported in Ref. l58l . 
Moreover, by fitting the size dependence of the numeri- 
cal results to the scaling law 1/L n , Franceschetti et a/.— 
reached n = 1.93 for InP nanocrystals and n = 1.97 for 
CdSe nanocrystals, whereas in Ref. [H, n = 2.51 was ob- 
tained for Si QDs. This discrepancy to the 1/L 3 law was 



8 



attributed to the presence of the LR component of the 
e-h exchange interaction. Obvious confusion is seen and 
hence a reexamination is necessary. 

We point out that not only the values of the exchange 
splittings, but also the scaling laws, depend on the dot 
size. So the investigation of size scaling of the e-h ex- 
change interaction is carried out in two limits: the strong 
and weak confinement limits. Moreover, 3D and 2D scal- 
ings are performed. For the 3D scaling, the dot size is 
varied in all three dimensions; whereas for the 2D scaling, 
the dot height is fixed and only the lateral size is varied. 
The characteristic size of the variation is denoted as L 
and the dot height l z in the 2D scaling is fixed and as- 
sumed to be much smaller than the exciton Bohr radius. 

Since the leading terms of the e-h exchange interaction 
originate from the heavy-hole-exciton basis in cubic III- 
V semiconductor QDs (e.g., GaAs, InAs and GaN QDs 
investigated in this paper, shown in Sec. Ill), we focus 
on these terms in Eqs. (1231) - (f3"0)) . i.e., terms without S*j 
and SL. Then the first terms in Eqs. (|2"4")) - (j2"u]) scale as 



dqUq( 



iq-(ri-r 



'QiQj 



did, 



4ire K\r[ - r 2 1 L 3 ' 



For the strong confinement limit where L is much smaller 
than the exciton Bohr radius, the exciton envelop func- 
tion F cv (ri , r 2 ) scales as oc -p in the 3D scaling and oc -p 
in the 2D scaling due to the normalization condition. 37 
In the calculation of matrix elements in Eq. @, we reach 
the scaling laws for the LR exchange terms 



111 



1 



TaTaTI oc Th (3D) 



L 3 L 3 L 3 
111 



L 3 
1 



w Vfp k f (2D) 

and for the short-range exchange terms 



1 1 



1 



(33) 
(34) 

(35) 
(36) 



These are the same as those in Ref. 1371 . 

However, for the weak confinement limit where L is 
much larger than the exciton Bohr radius, the relative 
motion of electron and hole is not sensitive to the value 
of L. Only the motion of the center-of-mass of the exciton 
is affected by the confinement. Therefore the scaling of 
the exciton envelope function changes to F c „(ri, r 2 ) oc -p 
in the 3D scaling and oc j- in the 2D scaling. Here we 
have assumed that l z is still smaller than the exciton Bohr 
radius in the weak confinement limit under consideration. 
Now terms of the LR exchange interaction scale as 



L 3 L 3 1 1 1 oc 1 

L L Z^Z^Z 



l2l2 T 3 ll*l 



ill 



l 



(3D), 
(2D), 



(37) 
(38) 



The short-range exchange terms scale as 

oil 1 , , ,11 1 , s 



(39) 



For genuine situation between these two limits, terms 
of the LR exchange interaction scale in the range of 
-3- ~ -p- in both 3D and 2D scalings while those of the 
short-range exchange interaction scale in the range of 

~ -p in the 3D scaling and -p ~ -p in the 2D scaling, 
respectively. In this way, we are able to explain the con- 
fusion discussed at the beginning of this subsection: the 
scaling of the e-h exchange interaction was investigated 
in different regimes of the confinement strength. We will 
further check these scaling rules in the next section in 
the analysis of the numerical results. 



C. Exciton spin relaxation 

The spin relaxation between the lowest two linear po- 
larized exciton states \X) and \Y) induced by the asym- 
metry of QD determines the dynamics of the optical lin- 
ear polarization decay. The \X) and \Y) states are de- 
fined as \X) = (|+1) + |-1)) A/2 and \Y) = -i(|+l) - 
\ — l))/\/2 where |±1) denotes the optically active exciton 
states with total angular momentum in the z-direction 
J z = ±1f2£ In GaN QDs, the effect of surface roughness 
can be suppressed experimentally to such an extent that 
it can be ignored. Therefore, the exciton spin relaxation 
is mainly assisted by electron- and hole-phonon interac- 
tions induced by deformation potential and piezoelectric 
field. From the Fermi golden rule, the phonon-assisted 
relaxation rate from \i) to |/) can be calculated by 



2vr 



- £ |M qA | 2 |</|x|i)| 2 [n qA <5( e/ - Ei - ft qA ) 



q,\ 



+ {nq\ + l)S(ef - €i + Hqx)], 



(40) 



in which M qA and \ 



gMj i-2 come f rom the 



electron- and hole-phonon interaction. n qA is the Bose 
distribution of phonon with mode A and wave vector q. 
In our calculation, the temperature is fixed at K, so the 
phonon absorption process is absent. 

In this paper, we take into account the electron- 
and hole-acoustic-phonon scattering due to the defor- 
mation potential with |M qs /| 2 = fiE. 2 q/(2pv s i), and 
due to the piezoelectric coupling with |A/ qs ;| 2 = 
288hir 2 e 2 el 4 (q x q y q z ) 2 /(K 2 pv s iq 7 ) for the longitudinal 
mode and E i=lj2 |M qptj | 2 = 32/i7r 2 e 2 e 2 4 [ g 2 g 2 + q 2 y q 2 z + 

lllx ~ (3<lxqyqz) 2 / q 2 }/ (k 2 pvstq 5 ) for the two transverse 
modest - — Here S , p, en and K stand for the acoustic 
deformation potential, the volume density of the mate- 
rial, the piezoelectric coupling constant and the static 
dielectric constant, respectively. v s i (v s t) is the longitu- 
dinal (transverse) sound velocity. Their values are given 
in Table I. 
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D. Numerical scheme 

In the computation we employ a disk-like QD model 
to simulate the real QDs, with infinite square-well po- 
tential in the z-direction (i.e., the [001] direction) and 
anisotropic coaxial harmonic-oscillator potential as the 
in-plane confinement! 83 ! 100 

t/ m ( 0, < z < l z . s 

V)W = (oo, others ' ^ 

V e (h)(x,y) = ^rn* em {ul e{h) x 2 + u 2 ye{h) y% (42) 

where u) xe{h) = ^*^r and ^(h) = T^f^J wi *h 
m *e(h)\\ denoting the effective mass of the electron (hole) 



in the plane. l x and l v are the characteristic lengths 
of the harmonic-oscillator potentials along the x- and y- 
directions and correspond to the major and/or minor di- 
ameters of the elliptic QD in the plane. l z corresponds 
to the dot height. In our model the single electron and 
hole experience different in-plane potentials but share the 
same confinement length. We adjust the relative magni- 
tudes of l x and l y to control the anisotropy of the QD 
and the magnitudes of l x , l y and l z to vary the strength 
of the confinement. 



TABLE I: Material parameters used in the calculation (from Ref. [73 unless otherwise specified). 





E a (eV) 


E p (eV) 


ml /mo 


m* o /m 


7i 


72 


K 


S(eV) 


GaN 


3.299 


25.0 


0.15 


0.29 


2.67 


0.925 6 


9.7 a 


8.3 d 


GaAs 


1.519 


28.8 


0.0665 


0.172 


6.85 


2.5 C 


12.53 a 


8.5 b 


InAs 


0.414 


21.5 


0.023 


0.14 


20.4 


8.7 e 


15.15 a 


5.8 b 




A SO (eV) p (10 3 kg/m 3 ) 


e M (10 8 V/m) 


AE LT (AteV) AE S r (fieV) 


v st (10 3 m/s) 


Vst (10 3 m/s) 




GaN 


0.017 


6.095 d 


43 d 






6.56 d 


2.68 d 




GaAs 


0.341 


5.31 b 


14.1 b 


80 b 


20 b 


5.29 b 


2.48 b 




InAs 


0.38 


5.9 b 


3.5 b 




0.3 C 


4.28 b 


1.83 b 





a Ref . Us, b Ref. 97, c Ref. .66, d Ref . M, 
° Obtained from Ref. [73] by 5(72 + 73) in the spherical approximation J'' 73 = 72 in this paper. 



The eigen equation for the envelope function 
F m n (1*1, T2) [Eq. @] is solved by the exact diagonal- 
ization method. The total Hamiltonian [Eq. (jAl[) ] is 
separated into two parts: H eh = H + H' where H 
[Eq. (|C24|) ] is diagonal in the real space and easy to be 
solved, and H' is the remaining parts of H eh which in- 
clude the e-h exchange interaction [Eqs. (|2"5)l - (I3TJ1) ] and 
H' given in Eq. (|C15[) . The eigenfunctions of Hq are 
taken as the basis functions, and the total Hamiltonian 
is diagonalized in the Hilbert space spanned by them. 
Detailed procedures are laid out in Appendix C. 

We stress that the direct Coulomb interaction is 
included in Hq, together with the confinement, and 
solved exactly. As pointed out in the Introduction, 
the direct Coulomb interaction is comparable or even 
stronger than the lateral confinement of the QDs un- 
der investigation*^ So it is not appropriate to con- 
struct the basis of exciton envelope functions by the 
product of single-particle wave functions of electron and 
hol o 69 ' 70 and treat the direct Coulomb interaction per- 
turbatively. As we will show below, the strength of the 



direct Coulomb interaction strongly affects the calculated 
fine structure splittings. Moreover, the big exciton bind- 
ing energy obtained^ indicates that the exciton wave 
functions, as well as the the spectrum, are markedly mod- 
ulated by the direct Coulomb interaction. Therefore, the 
unperturbative treatment of the direct Coulomb interac- 
tion is crucial to fully take into account its effect on the 
exciton fine structures. 



III. EXCITON FINE STRUCTURE IN QDS 

In this section, we investigate the exciton fine structure 
in single GaAs, InAs and GaN QDs. The fourfold degen- 
erate exciton ground states are split into a J z = ±2 dark 
doublet and two bright states when the e-h exchange in- 
teraction is taken into account. In a circular QD, the two 
bright states are degenerate with the ^-component of to- 
tal angular momentum J z = ±1. For an anisotropic QD, 
the e-h exchange interaction couples these two states to- 
gether, forming the so-called \X) and \Y) exciton states. 
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These exciton fine structures are schematically shown in 
Fig. [TJ The dark exciton doublet does not split since 
the e-h exchange interaction is absent for J z — ±2 states 
[Eqs. (g31) and ([27])]. 103 lt is noted that in our results the 
exciton ground state is always the dark doublet, which 
is consistent with the common understanding . 19 ' 38 ! 53 We 
restrict ourselves in discussing the fine structures split 
from the originally fourfold-degenerate exciton ground 
states when the e-h exchange interaction is introduced, 
unless otherwise specified. 



SR 



SR+LR 



Isotropic QD 4-fold 



Anisotropy QD 4-fold 



-<±i> 
-<±2> ■ 



<±l> 

BD exchange splitting 
<±2> 



/ K }l*>.|Y> 
x>,|Y>^"- ! — \ J 



= |X>,|Y> 

• <±2> 



Doublet splitting energy 
<±2> 



FIG. 1: The exciton fine structures in isotropic and 
anisotropic QDs are shown schematically for the cases with 
no e-h exchange interaction (left column), with only the SR 
exchange interaction (middle column) and with both the LR 
and SR exchange interactions (right column). Note that the 
doublet splitting energy with only the SR exchange interac- 
tion is very small and has been exaggerated in the figure. 



A. GaAs and InAs 

We first study the exciton fine structures in GaAs and 
InAs QDs. The doublet splitting energy and the BD 
exchange splitting are calculated and scaling analyses are 
performed. 



direct Coulomb interaction. When the direct Coulomb 
interaction is totally switched off, the doublet splitting- 
energy is only less than 10% of its original value where 
the direct Coulomb interaction is fully taken care of. 
The strong dependence of doublet splitting energy on 
the strength of direct Coulomb interaction is explained 
as follows. The direct Coulomb interaction attracts the 
electron and hole together and enhances the overlap of 
their wave functions. Hence according to Eqs. (f25j) and 
(J2U), with stronger direct Coulomb interaction, stronger 
e-h exchange interaction is obtained. 

Another important feature in Fig.[2]is that the doublet 
splitting energy strongly depends on the dot shape. The 
absolute value of the doublet splitting energy decreases 
with decreasing dot anisotropy, and tends to zero when 
the confining potential approaches isotropic. It is seen 
from the curve with 7/ = 1 that the doublet splitting 
energy varies from to about 250 /ieV when l y is varied 
from 20 nm(=/ x ) to 10 nm. 

Our results are much larger than those reported very 
recently by Kadantsev and Hawrylak,— where they cal- 
culated the doublet splitting energy with exciton wave 
function constructed by the product of single-particle 
ground-state wave functions of electron and hole states. 
So the effect of the direct Coulomb interaction was to- 
tally absent and the doublet splitting was hence underes- 
timated. Detailed comparison with the results in Ref. [6^ 
is given in Appendix D. 

Our results are in good agreement with the existing 
experimental results! 21 ' 29 ' 104 For example, the doublet 
splitting energy measured by Gammon et al. 21 in sin- 
gle GaAs QD lies in the range 20~50 /xeV, which cor- 
responds approximately to l y = 16~18 nm, l x — 20 nm 
and l z — 3 nm in our model as shown in Fig. [5J Good 
agreement is also reached with former theoretical work by 
Takagahara— based on the variational method with the 
direct Coulomb interaction included, which also showed 
good agreement with the same experiment J^ 5 - 



1. Doublet splitting 



In order to investigate the influence of the direct 
Coulomb interaction on the doublet splitting energy, we 
introduce a dimensionless parameter r\ in front of the di- 
rect Coulomb interaction U ch in Eq. (jAll) . By varying 
rj from to 1, the direct Coulomb interaction is varied. 
In the calculation, the major/minor diameter along the 
x-direction l x is fixed at 20 nm and the minor/major 
diameter along the y-direction l y is varied from 10 to 
30 nm. The disk-like dot height l z is fixed at 3 nm. In 
our model, the ground state of bright exciton is found to 
polarize along the axis of the weaker confinement, i.e., 
| A) state when l x > l y and \Y) state when l x < l y . 
Here the \X) and \Y) states are defined by their main 
components. The doublet splitting energy is defined as 
Ey — Ex here and hereafter, with Ex and Ey repre- 
senting the eigenenergies of \X) and \Y) exciton states, 
respectively. From Fig. [21 one observes that the doublet 
splitting is markedly reduced with the decrease of the 



TABLE II: Size parameters used in the calculation of doublet 
splitting energy in single InAs QD for Fig. [3jb) . 





1 


2 


3 


4 


5 


6 


7 


lx (nm) 


21 


19 


17 


15 


14 


13 


12 


l z (nm) 


5.5 


5.4 


5.3 


5.2 


5.1 


5 


4.9 




8 


9 


10 


11 


12 


13 




lx (nm) 


11.5 


11 


10.5 


10 


9.5 


9 




l z (nm) 


4.8 


4.75 


4.7 


4.55 


4.5 


4.4 





Similar features of the shape dependence of doublet 
splitting energy are shown in Fig. [31(a) for InAs QDs 
with the direct Coulomb interaction fully included. In 
addition, we plot the doublet splitting energy as a func- 
tion of the exciton recombination energy in Fig. [3]Jb) in 
order to compare with the experimental data by Seguin 
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FIG. 2: (Color online) The doublet splitting energies in single 
GaAs QD as function of the dot major/minor diameter l y 
with fixed minor/major diameter l x = 20 nm and dot height 
l z = 3 nm. Different values of r\ are chosen to modulate the 
strength of the direct Coulomb interaction. 



et alQ. The exciton recombination energy is defined 
as Ex + Ey + Eg, with E g standing for the band gap 
in cubic InAs. The experimental data are taken from 
Ref. [3ll In the calculation, we fix the major/minor di- 
ameter l y = 10 nm and vary the minor/major diameter 
l x in the range of 8.5 ~ 20 nm. The dot height l z is 
varied in the range of 4.3 ~ 5.5 nm. The exciton re- 
combination energy is mainly modulated by the strong 
confinement along the z-direction while the doublet split- 
ting energy is mainly determined by the ratio of l x : l y . 
So the theoretical points are obtained by properly choos- 
ing the values of (l x , l z ) in pairs. From left to right in 
Fig. [3[b), the explicit values of (l x , l z ) for the theoreti- 
cal points are given in order in Table II. It is seen from 
the figure that our theoretical results are in good agree- 
ment with the experimental data. The doublet splitting 
energy decreases with increasing exciton recombination 
energy and intersects through zero. Generally speaking, 
we employ larger values of l x together with larger l z , and 
hence smaller recombination energy. The zero point of 
doublet splitting energy is reached at l x = 10 nm and 
l z = 4.55 nm. In this way, the trend of the variation 
of the doublet splitting energy with the exciton recom- 
bination energy is well explained as a result of the dot 
geometry*^ 

The relative importance of the LR and SR exchange 
interactions with respect to the exciton fine structure 
splittings in QDs is an open question. Early work by 
Efros et al. included only the SR exchange interac- 
tion to investigate the band-edge excitons in spherical 
QDs. Takagahara2£ assigned the origin of exciton dou- 
blet structure to the LR exchange interaction and Glazov 
et alQ took into account only the LR exchange interac- 
tion in their work to investigate exciton fine structure 




Theory 
Experiment 



(b) InAs 



1.05 1.1 1.15 1.2 1.25 1.3 1.35 
Recombination energy (eV) 



FIG. 3: (Color online) The doublet splitting energy in single 
InAs QD (a) as a function of the dot major/minor diameter l v 
with fixed minor/major diameter l x = 20 nm and dot height 
lz — 3 nm; (b) as a function of the exciton recombination 
energy. The experimental points are taken from Ref. |3l| and 
the theoretical points are calculated with the size parameters 
listed in Table II. The dotted line is plotted for the guide of 
eye. 



in an anisotropic QD. However, Tsitsishvili et a?j 66 i 72 
and Horodyska et al. in their very recent works! took 
into account only the SR exchange interaction for exci- 
tons in anisotropic and spherical QDs, respectively. Here 
we reexamine the relative importance of the LR and 
SR exchange interactions to the doublet splitting energy. 
Those to the BD exchange splitting are to be discussed 
in the following corresponding parts. 

We calculate the doublet splitting energies by includ- 
ing the LR, SR and both LR and SR exchange inter- 
actions. In Fig. [H we plot the doublet splittings as a 
function of the major/minor diameter l y with the mi- 
nor/major diameter l x fixed in GaAs and InAs QDs. It 
is seen that the curves obtained by including the LR ex- 
change interaction only and by including both the LR and 
SR exchange interactions almost match each other. The 
doublet splitting energies from the SR exchange interac- 
tion are always less than 1 /ieV, more than two orders 
of magnitude smaller than the splittings caused by the 
LR exchange interaction. So the LR exchange interaction 
is dominant in determining the doublet splitting energy 
when all heavy-hole, light-hole and split-off bands are 
taken into account. This is consistent with Ref. 1371 where 
only the heavy- and light-hole bands are taken account 
of. ' 

It is also noted that the doublet splitting energies from 
the SR exchange interaction are in opposite sign to those 
from the LR exchange interaction. This means that the 
relative positions of \X) and \Y) exciton states are re- 
versed when only the SR exchange interaction is included. 

We further investigate the size dependence of the dou- 
blet splittings of GaAs and InAs QDs and perform the 
3D and 2D size-scaling analysis. Detailed results are plot- 
ted in Fig. [5] with the solid curves fit to the power law 
oc 1/L n . The results are consistent with the physical 
intuition. The doublet splitting energy increases with 
the decrease of the dot size since for smaller dot size, 
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TABLE III: Power indices of the 2D and 3D scalings of the 
BD exchange splitting contributed by the LR, SR and both 
LR and SR exchange interactions in GaAs and InAs QDs. 





3D 


2D 




LR 


SR 


Both 


LR 


SR 


Both 


GaAs 


1.57 


1.52 


1.49 


1.14 


0.15 


0.22 


InAs 


1.95 


2.00 


1.95 


1.69 


0.80 


1.38 



2. BD exchange splitting 



FIG. 4: (Color online) The doublet splitting energies from 
the LR, SR and both LR and SR exchange interactions, are 
plotted as a function of the dot major/minor diameter l y , 
with fixed minor/major diameter l x = 20 nm and dot height 
h = 3 nm for (a) GaAs QD and (b) InAs QD. 



the overlap of the wave functions of electron and hole 
is enhanced and hence larger matrix elements of e-h ex- 
change interaction are obtained [Eqs. (1251) and (|29|) and 
their Hermitian conjugates]. We obtain TiGaAs = 1-45 
and ni n As = 2.14 in the 3D scaling and iiQaAs = 1-04 and 
^inAs = 1-6 in the 2D scaling. All power indices obtained 
lie in the range of 1^3 and are consistent with the scal- 
ing rules for the LR exchange interaction established in 
Sec. II. The fact that the power indices of InAs QDs are 
larger than those of GaAs QDs is understood as follows. 
The exciton Bohr radius is 14.9 nm in bulk GaAs and 
51.6 nm in bulk InAs. As compared to the characteristic 
size L of the QDs in the 3D and 2D scalings, for GaAs 
QDs, ag^ r s is comparable to or smaller than L. So GaAs 
QDs are closer to the weak confinement limit. For InAs 
QDs, a^^. is larger than ago^f, so InAs QD is closer to 
the strong confinement limit compared to GaAs QD. 
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FIG. 5: (Color online) The doublet splitting energies in GaAs 
and InAs QDs. (a) For the 3D scaling, with l x : l y : l z = 6 : 4 : 
1. (b) For the 2D scaling, with l x : l y — 4 : 3 and l z = 3 nm. 
The solid curves are fit to the power law C/L n and the index 
values n are shown next to the curves. 



The band-edge exciton states in isotropic QDs were 
investigated in previous work o 65 i 67 by including only the 
SR exchange interaction. It is noted that even though 
the doublet splitting energy, which is mainly from the LR 
exchange interaction, becomes zero in circular QDs, the 
LR exchange interaction still contributes to the splitting 
between the bright and dark exciton states. For circular 
GaAs QD with l z = 3 nm and l x = l y = 20 nm, the 
BD exchange splitting is 1.76 meV and the LR exchange 
interaction contributes 0.17 meV to it. So in GaAs QDs, 
the SR exchange interaction is dominant in determining 
the BD exchange splitting. It is a different case in InAs 
QD of the same size, where the LR exchange interaction 
contributes 0.31 meV out of the total BD exchange split- 
ting 0.49 meV. This difference is due to the fact that 
the singlet-triplet splitting parameter AE'sr is 20 fieV in 
GaAs compared 0.3 /LteV in InAs. As a result, the SR 
exchange interaction is rather weak in InAs. The calcu- 
lated BD exchange splitting of 500 /ieV in InAs QD with 
l z = 6 nm and l x = ly = 15 nm is in good agreement 
with the experimental data^ So we emphasize that the 
LR exchange interaction is important to understand the 
experimental results not only for the doublet splitting 
energy, but also for the BD exchange splitting in QDs. 

The size dependence of the BD exchange splitting is 
investigated in circular QDs for the purpose of eliminat- 
ing the effect of doublet splitting. In the calculation, 
we take into account the LR, SR and both LR and SR 
exchange interactions. The results are shown in Fig. [5] 
The solid curves are fit to the power law cx l/L n , and 
the obtained power indices of the size scaling of BD ex- 
change splitting are listed in Table III. The size scaling 
laws of the BD exchange splitting by including only the 
LR exchange interaction are close to those of the dou- 
blet splitting energy in both the 3D and 2D scalings in 
GaAs and InAs QDs (see Fig. M and Table III). This is 
because the diagonal and off-diagonal matrix elements of 
the LR exchange interaction in Eq. (|23[) are in similar 
forms of the r dependence. The power indices of the BD 
exchange splitting from the SR exchange interaction lie 
reasonably in the range 1~3 in the 3D scaling and in the 
range 0~2 in the 2D scaling, and are consistent with the 
scaling rules established in Sec. II. 

An important feature of the relative energy positions of 
the dark and bright exciton states can be deduced from 
the above scalings. As in the 2D scalings, the power 
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FIG. 6: (Color online) The BD exchange splittings in circular 
QDs by taking account of the LR, SR and both LR and SR 
exchange interactions, (a) For the 3D scaling in GaAs QD, 
with l x : l y : l g = 6 : 6 : 1. (b) For the 2D scaling in GaAs 
QD with l x — l y and l z = 3 nm. (c) For the 3D scaling in 
InAs QD, with l x : l y : l z = 6 : 6 : 1. (d) For the 2D scaling 
in InAs QD, with l x = l y and l z = 3 nm. The solid curves are 
fit to the power law C / L n with n listed in Table III. 

indices of the BD exchange splitting in isotropic GaAs 
and InAs QDs (also in GaN QD, see in the next subsec- 
tion) are smaller than 2 in the size range under inves- 
tigation. Meanwhile the level spacing between ground 
dark level and first excited dark level, which is induced 
by the lateral confinement, scales approximately as cx -p- 
[Eq. (|C21[) ].i££ So as the dot diameter increases, the first 
excited dark exciton level decreases more rapidly than 
the ground bright exciton level, and a crossing between 
these two levels may occur. For a typical case, in Fig. [7] 
we plot the eigenenergies of the ground dark exciton and 
bright exciton levels and the first excited dark exciton 
level as a function of the dot diameter l x (=l y ) in GaAs 
QD. The dot height l z is fixed at 3 nm. When the dot di- 
ameter l x is varied from 10 to 30 nm, a crossing between 
the ground bright doublet and the first excited dark quar- 
tet (since l x = l y , the first excited dark exciton level is 
fourfold-degenerate) is observed around l x = 23 nm. This 
crossing can also be obtained in anisotropic QDs with the 
underlying physics. 

This size-dependent bright-dark exciton level cross- 
ing provides a unique way of tuning the bright exciton 
states in resonance with the dark exciton states, which 
is meaningful to recent research on the optical nuclear 
spin pumping with the help of the hyperfine-interaction- 
mediated spin-flip transitions between the bright and 
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FIG. 7: (Color online) The energy positions of the lowest dark 
and bright levels and the first excited dark level as function 
of l x , with l x — l y and l z — 3 nm. The zero point of energy is 
set at the lowest energy level at l x = 30 nm. 



B. GaN 

We then investigate the properties of the LR and SR 
exchange interactions in cubic GaN QD where theoretical 
works are absent. Compared to Ini-^Ga^ As-based struc- 
tures, nanostructures based on GaN are less investigated 
both experimentally and theoretically. Much attention 
has been attracted for their unique properties, e.g., the 
wide band gap (^3.3 eV), which represents great poten- 
tial for applications in electronics and photonics at tem- 
perature much higher than the liquid-helium or liquid- 
nitrogen cryogenic temperatures! 39 ' 1 10 ' 1 1 1 

1. Doublet splitting 

Due to the absence of experimental value of the singlet- 
triplet splitting A_Esr in cubic GaN, the doublet splitting 
energy in GaN QDs is calculated by temporarily setting 
A^sr = 20 fieVm In Fig. Ufa), we plot the doublet 
splitting energies calculated by including the LR, SR and 
both LR and SR exchange interactions as a function of 
the dot minor diameter l y . Similar to GaAs and InAs 
QDs, the SR exchange interaction is irrelevant when con- 
sidering the doublet splitting energy. Since the strength 
of the SR exchange interaction is proportional to the 
value of AEsr, we assert that ASsr in GaN can not 
be so large as to make the SR exchange interaction com- 
parable or even exceed the LR exchange interaction be- 
cause otherwise it will make other results, e.g., the BD 
exchange splitting, unreasonable. So in Fig.[5fb), we plot 
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the doublet splitting energy varying with the dot shape, 
calculated without the SR exchange interaction. Size pa- 
rameters are chosen according to the experiment.— One 
observes in the figure that with large anisotropy, the dou- 
blet splitting energy reaches 100s of /zeV. The large dou- 
blet splitting energy obtained is key to understand the 
experiment by Lagarde et al^ where the conversion from 
exciton linearly polarized states \X) and \Y) to the circu- 
larly polarized ones | ± 1) was not observed for magnetic 
field up to 4T. 
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FIG. 8: (Color online) The doublet splitting energies in single 
GaN QD as a function of the dot major/minor diameter l y 
with fixed minor/major diameter l x = 12 nm and dot height 
l z — 2 nm. (a) The doublet splitting energies from the LR, 
SR exchange interactions and both together; (b) the doublet 
splitting energy from the LR exchange interaction. 

The size dependence of doublet splitting energy is also 
investigated. The obtained results can be well under- 
stood in a straightforward way as those for GaAs and 
InAs QDs discussed above. So we only plot our results 
in Fig. [9] without more discussions. The power indices 
obtained from the 2D and 3D size scalings are shown 
next to the curves. The scaling laws of doublet splitting 
energies found in GaN QDs are consistent with the rules 
established in Sec. II. 

We further calculate the doublet splitting energy in 
small cubic GaN QDs by pushing our model to its 
extreme.*^ In Fig. [TUl we plot the doublet splitting en- 
ergies as function of dot size. It is seen that the doublet 
splitting energies reach several cV when the QDs become 
extremely small. This is strongly supported by recent 
experiment on wurtzite GaN QDs where doublet split- 
ting energies in the range of 2^7 meV were reported. 39 
Further experiments on cubic GaN QDs are expected. 



2. BD exchange splitting 

Although the SR exchange interaction is negligible con- 
cerning the doublet splitting in GaN QDs, its contribu- 
tion to the BD exchange splitting is still comparable to 
that of the LR exchange interaction. In order to investi- 
gate the properties of BD exchange splitting in GaN QDs, 
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FIG. 9: (Color online) The doublet splitting energies in GaN 
QDs. 1) For 3D scaling, with l x : l y : l z = 6 : 4 : 1. l z is varied 
from 1.5 to 3.5 nm. 2) For 2D scaling, with l x : l y = 4 : 3 and 
l z = 2 nm. The solid curves are fit to the power law C / L n 
with the index values n shown next to the curves. Note the 
scale of the 2D scaling is on top of the frame. 
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FIG. 10: (Color online) The doublet splitting energies in small 
cubic GaN QDs as function of dot size. Note the scale of one 
of the curves is on top of the frame. 



the singlet-triplet splitting AEsn is again set at 20 /ieV 
as a parameter Due to the fact that in circular QDs, 
the BD exchange splitting calculated by including both 
the LR and SR exchange interactions is approximately 
the summation of those obtained by including each sep- 
arately and the contribution of the SR exchange interac- 
tion is proportional to the value of AE'sr,, the genuine 
value of AiJgp, in cubic GaN can be extracted by com- 
paring our theoretical results of BD exchange splitting 
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TABLE IV: Power indices of the 2D and 3D scaling of the 
BD exchange splitting contributed by the LR, SR and both 
LR and SR exchange interactions in GaN QDs. 
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GaN 
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1.26 


0.22 
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with future experimental data. 

In Fig. [TTJ we plot the BD exchange splitting as func- 
tion of the dot size. The obtained power indices of size 
scaling are listed in Table IV. From Fig. [LTTbL one ob- 
serves a crossing of BD exchange splittings obtained from 
the LR exchange interaction and from the SR exchange 
interaction, indicating the exchange of the relative impor- 
tance of the LR and the SR exchange interactions. This 
is due to the different scaling rules of the LR and SR 
exchange interactions in the 2D scaling [Eqs. (|34|). (|38|) 
and (j3l)]) . (j3"9|) ]. As the the dot diameter increases, the 
BD exchange splitting from the LR exchange interaction 
decreases more rapidly than that from the SR exchange 
interaction (see in Table IV for detailed values of power 
indices) and becomes smaller than it for l x > 11.4 nm 
in GaN QD. In fact, one may also expect a crossing in 
Fig. HJd) in the region of l x > 30 nm in InAs QD. 
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FIG. 11: (Color online) The BD exchange splittings in GaN 
QDs by taking account of both the LR and SR exchange 
interactions, and separately: (a) for the 3D scaling, with 
lx '■ l y : l y — 6 : 6 : 1 and (b) for the 2D scaling, with 
lx — ly an d l z — 2 nm. The solid and dashed curves are fit to 
the power law C/L n and the index values n are shown next 
to the curves. 



C. Discussions on relative importance of the 
heavy-hole, light-hole and split-off bands 

We now turn to address the relative importance of 
the heavy-hole, light-hole and split-off bands to the ex- 
citon fine structure. Theories of exciton fine structure 
in QDs were presented by taking into account only the 
heavy-hole band i 62 ' 115 It was pointed out by Takagahara 
that it is important to take account of the mixing of 
heavy-hole and light-hole bands to explain the experi- 



mental phenomenal But this argument was not explic- 
itly proved in that paper. In Fig. [T^l we plot the en- 
ergies of the lowest subbands of decoupled heavy, light 
and split-off holes in GaAs, InAs and GaN QWs. The 

^2 2 

valence band coupling is neglected so that E z — , 
with m* representing the effective mass of the heavy, light 
or split-off hole in the z-direction. The relative positions 
of these three subbands vary with l z because of different 
hole effective masses in the z-directionjZ2 As shown in 
the figure, in the strong confinement regime, the split- 
off subband is closer to the heavy-hole subband than the 
light-hole one for all three materials. Therefore, if the 
light-hole band is considered, the split-off band should 
also be included. 

In order to investigate the relative importance of the 
heavy-hole, light-hole and split-off bands to the e-h ex- 
change interaction, the doublet splitting energy and the 
BD exchange splitting are calculated by taking into ac- 
count: (i) the heavy-hole band only; (ii) heavy-hole and 
light-hole bands; and (iii) all the three valence bands, 
separately. In Fig.[T31 we plot the doublet splittings orig- 
inated from the three valence bands separately. The con- 
tribution of the heavy-hole band is calculated by includ- 
ing the terms of e-h exchange interaction derived from 
the heavy-hole band only as in Case (i). The contribu- 
tion from the light-hole (split-off) band is obtained by 
subtracting the splitting calculated from Case (i) [Case 
(ii)] from the one from Case (ii) [Case (iii)]. It is noted 
that the doublet splitting energy is decreased by fur- 
ther including the light-hole and split-off bands. There- 
fore in Fig. [131 the contributions from the light-hole and 
split-off bands are multiplied by —1. Moreover, for GaN 
QDs, the structure splittings are calculated by setting 
A-Esr = 20 neV^m As shown in the figure, the contri- 
bution of the heavy-hole band is much larger than those 
from the other two bands. The doublet splitting energy 
is slightly changed by further including the light-hole and 
split-off bands as in Case (ii) and (iii), with the split-off 
band contributing most of the change. Under all size 
parameters adopted, the doublet splitting energy is only 
changed by less than 2 % in both GaAs and GaN QDs, 
and less than 7.2 % in InAs QDs. One can hence conclude 
that the terms derived from the heavy-hole band domi- 
nate the off-diagonal matrix elements of the LR exchange 
interaction [Eq. (1231) and its Hermitian conjugate], i.e., 
\Sij\ <C 1. And they are also much larger than the off- 
diagonal matrix elements of the SR exchange interaction 
[i.e., Eq. (|2T)]) and its Hermitian conjugate]. 

Similarly, we find that the contribution of light-hole 
and split-off bands to the BD exchange splitting is also 
small, in the order of 1 /xeV, which is much smaller than 
100s of fj,eV from the heavy- hole band for all three ma- 
terials under investigation. This can be understood in 
the same way: according to Eqs. (l2Ul) - (j3"0)l . for the LR 
exchange interaction, terms without Sij and Si* should 
also dominate the diagonal matrix elements of the LR ex- 
change interaction which contribute to the BD exchange 
splitting. Moreover, for the SR exchange interaction, the 
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diagonal matrix elements [Eqs. (|28|) and (|30|) ] are not 
affected by the inclusion of the light-hole and split-off 
bands up to the order under consideration. 

In short, both the light-hole and split-off bands are 
negligible when the band-edge exciton fine structure 
is investigated in cubic III-V semiconductor QDs with 
strong confinement along the [001] direction. This fur- 
ther demonstrates the feasibility of treating the light-hole 
and split-off bands perturbatively through the Lowdin 
partition. However, as discussed after Eqs. (|21l) and (|22p. 
the confinement-induced valence band mixing is crucial 
in understanding some exciton properties, e.g., the ob- 
servability of the dark excito n 28 i 38 ' 87 ' 88 and the degree of 
the linear polarization of QD emission.*!^ 




\ (c) GaN 


HH 




LH 




SO 



1.5 2 2.5 3 3.5 4 4.5 5 2 3 4 5 6 7 8 9 10 1.5 2.5 3.5 4.5 5.5 

l z (nm) 



FIG. 12: (Color online) Eigenenergies of decoupled heavy-hole 
(HH), light-hole (LH) and split-off- hole (SO) ground states in 
QW. Infinite square well potential is employed and the well 
width is denoted as l z . (a) GaAs, (b) InAs, (c) GaN. 



from the Fermi golden rule with the exciton eigenfunc- 
tions and eigenenergies obtained from the exact diagonal- 
ization method. 116 Since the ground state of the bright 
exciton is polarized along the major axis of the potential 
ellipse, the relaxation rate is calculated from the upper 
state of the doublet to the lower one, i.e., |y) — > \X) for 
l x > ly and \X) -> \Y) for l x < l y . 

In the calculation, we fix the major/minor diameter 
l x = 12 nm and the dot height l z = 2 nniiS The shape de- 
pendence of the exciton spin relaxation rate is studied by 
varying the minor/major diameter l y and the results are 
plotted in Fig. Q31 One observes from the figure that the 
relaxation rate between the lowest \X) and \Y) exciton 
states shows strong dependence on the dot anisotropy. In 
the regime of large anisotropy, the relaxation rate reaches 
10 4 /xs _1 , which in turn gives the exciton spin relaxation 
time in the order of 10 ps. In the vicinity of l y = l x , 
the relaxation rate decreases quickly with decreasing dot 
anisotropy. When the QD approaches the circular shape, 
the rate of change is even larger. The relaxation rate 
decreases down to less than 100 s _1 and tends to zero 
when l y approaches l x . This calculated long exciton spin 
relaxation time, especially in the range of l y 11~13 nm, 
is supported by the latest experiment^ 

It is pointed out that the behavior of spin relaxation 
rate in the range of l y — 7 ~16 nm mostly results from 
the anisotropy dependence of doublet splitting energy 
shown in Fig. HJb). From Eq. (|40l) . one has 
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FIG. 13: (Color online) Contributions of the heavy-hole (HH), 
light-hole (LH) and split-off (SO) bands to the doublet split- 
ting energy, calculated as a function of the dot minor diam- 
eter l y , with fixed dot major diameter l x = 20 nm and dot 
height l z — 3 nm for (a) GaAs QD, (b) InAs QD; and with 
l x — 12 nm and l z = 2 nm for (c) GaN QD. 



IV. EXCITON SPIN RELAXATION IN GAN 
QDS 

In this section, we study the exciton spin relaxation 
in single GaN QDs. The relaxation rate is calculated 



When the wave vector q\ is not too large, the value of 
the relaxation rate is mainly modulated by the factor 
q^|M q ^| 2 , which increases with increasing q\ for all three 
channels under consideration. The wave vector of the 
acoustic phonon is given by q\ — j^-, with v\ repre- 
senting the sound velocity and AE denoting the phonon 
energy, which is equal to the difference of eigenenergies 
between the two bright exciton states considered, i.e., the 
doublet splitting energy. One observes from Fig. Mjo) 
that qualitatively, the doublet splitting energy is pro- 
portional to the dot anisotropy. So large dot anisotropy 
indicates large phonon wave vector, which in turn results 
in the anisotropy dependence of relaxation rate as shown 
in Fig. [TH in the range of l y = 7 ~ 16 nm. 

With further decrease of l y down to less than 7 nm, the 
exciton spin relaxation rate reaches a maximum around 
l y = 6 nm, where the doublet splitting energy is about 
730 fieV which corresponds to the situation that the 
wavelength of the emissive phonon becomes comparable 
with twice of the lateral dot sizei 117 ' 118 On the other side, 
when l y increases over 16 nm, the increase of the relax- 
ation rate becomes slower and when l y reaches 24 nm, 
the relaxation rate also reaches a maximum. The under- 
lying physics is different from the previous one. Here, 
the maximum is attributed to the interplay of the dot 
anisotropy and the strength of the lateral confinement. 
On one hand, the increase of l y gives a more anisotropic 
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dot shape which tends to increase the doublet splitting 
energy; on the other hand, the confinement along the y- 
direction becomes weaker with larger l y , which tends to 
decrease the overlap of the electron and hole wave func- 
tions and reduce the doublet splitting energy according 
to Eq. (|25|). The competing leads to the maximum. 
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FIG. 15: (Color online) The relaxation rate limited by the 
electron/hole- longitudinal acoustic phonon scatterings due to 
the deformation potential (DA) and the piezoelectric coupling 
(LA), and by the electron/hole-transverse acoustic phonon 
scattering due to the piezoelectric coupling (TA). The relax- 
ation rates are plotted as function of the dot major/minor 
diameter l y with fixed minor/major diameter l x = 12 nm and 
dot height l z = 2 nm. 



FIG. 14: The dependence of relaxation rate between the low- 
est \X) and \Y) bright exciton states on the dot major/minor 
diameter l y with fixed minor/major diameter l x — 12 nm and 
dot height l z — 2 nm. 

In addition, the relative importance of the three chan- 
nels contributing to the exciton relaxation rate is inves- 
tigated. We take into account one channel at a time and 
the results are plotted in Fig. [15] It is seen that the 
relaxation rate limited by the electron/hole-longitudinal 
acoustic phonon scattering due to the piezoelectric field 
coupling is always 1 to 3 orders of magnitude larger than 
that due to the deformation potential, but is 2 to 3 or- 
ders of magnitude smaller than the one limited by the 
electron/hole-transverse acoustic phonon scattering due 
to the piezoelectric coupling. So the transverse acoustic 
phonon-emission process dominates the exciton spin re- 
laxation between the lowest \X) and \Y) bright exciton 
states in GaN QDs. 



V. SUMMARY 

In summary, we have established a general scheme to 
investigate the exciton fine structure and spin relaxation 
in cubic III-V semiconductor QDs. A 12 x 12 matrix 
representation of the exciton Hamiltonian corresponding 
to the LR and SR exchange interactions is derived by 
taking into account the conduction band Tg, the heavy- 
hole and light- hole bands Tg and the split-off band T^, 
where the split-off band has never been included explic- 
itly to investigate the exciton properties. In the case 
with strong confinement in one direction (the [001] di- 
rection in this paper), the Lowdin partitioning method 



is employed to take account of the confinement-induced 
band mixing and a four-band Hamiltonian of the e-h ex- 
change interaction is derived. We use these formulae in 
the study of the relative importance of the heavy-hole, 
light-hole and split-off bands to the exciton fine struc- 
ture. We find that the contribution of the split-off band 
is a little larger than that of the light-hole band, but both 
are negligible when considering the exciton fine structure 
in GaAs, InAs and GaN QDs. This behavior in GaN 
QDs is unexpected since a significant effect of the split- 
off band is expected due to the large band gap and the 
small spin-orbit splitting in cubic GaN. We attribute this 
to the confinement-induced subband splitting due to the 
different effective masses of the heavy, light and split-off 
holes in the z direction. 

In our approach, the direct Coulomb interaction is 
treated unperturbatively. We find that the strength of 
the direct Coulomb interaction strongly affects the dou- 
blet splitting energy (hence also the BD exchange split- 
ting). We also show that previous works in which the 
direct Coulomb interaction was treated perturbatively 
vastly underestimate the doublet splitting. We demon- 
strate that the exact inclusion of the direct Coulomb in- 
teraction is important for excitons in the weak confine- 
ment regime .^2^ 

We also discuss the size and shape dependences of the 
doublet splitting energy and the BD exchange splitting. 
Strong anisotropy dependence of the doublet splitting en- 
ergy is reported, which agrees with the former theoretical 
and experimental works on GaAs and InAs QDs i 31 i 37 The 
size dependences of the doublet splitting energy and the 
BD exchange splitting are investigated by performing the 
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size-scaling analysis. The behavior of the variation of the 
fine structure splittings with the dot size is well explained 
by the scaling rules established. The doublet splitting 
energy in cubic GaN QDs increases with the increase of 
dot anisotropy and/or the decrease of dot size, varying 
from to 100s of /j.eV, and reaching up to several meV 
for extremely small dot size and large dot anisotropy. 
Our results are well supported by recent experimental 
finding o 39 ' 78 but call for more experimental works. The 
still undetermined singlet-triplet splitting Ai?gR in cubic 
GaN can be fit from future experimental data on the BD 
exchange splitting, but the uncertainty does not affect 
the conclusions of this paper. 

We investigate the relative importance of the LR and 
SR exchange interactions to the exciton fine structure. 
The LR exchange interaction is identified as the origin of 
the exciton doublet structure, which is in agreement with 
that by Takagahara 37 where only the heavy- and light- 
hole bands were included. We show that the LR exchange 
interaction, which is absent in many previous works ; 65 ! 67 
contributes to the splitting between the bright and dark 
exciton states, even in circular QDs where the doublet 
splitting vanishes. The contribution of the LR exchange 
interaction to the BD exchange splitting is smaller than 
that of the SR one in GaAs QDs but is comparable in 
InAs and GaN QDs. Our calculations also demonstrate 
that the relative importance of the LR and SR exchange 
interactions to the BD exchange splitting can exchange 
with the variation of the dot lateral size. 

The exciton spin relaxation in cubic GaN QDs is also 
investigated. We find that the exciton spin relaxation 
rate strongly depends on the dot anisotropy (relaxation 
rate from 10 4 /is -1 down to less than 10~ 4 ^s _1 is re- 
ported). In the small anisotropy regime, the long exciton 
spin relaxation time obtained (longer than 100 ns) is in 
good agreement with recent experiment by Lagarde et 
al.Z& The electron/hole-transverse acoustic phonon scat- 
tering due the piezoelectric field is recognized as the dom- 
inant magnetism of the exciton spin relaxation. 

Finally, we address the possible extensions of our work: 
(a) Other than the shape anisotropy of QDs, strain 
anisotropy can also result in the doublet splitting, which 
is not included in the present investigation, (b) Within 
our model, the infinite square well potential is employed 
as the confinement along the z direction and only the 
lowest elect ron and hole subband is included. As pointed 
out in Ref. Ill4 in the case of extremely small QDs, one 
needs to switch to other model potential or approach such 
as the psudopotential approximatio n 3 ^ 1 19 1 1 20 to obtain 
more accurate results. Meanwhile, for QDs, the vertical 
height of which is not much smaller than its lateral size 
or the exciton Bohr radius, the multi-subband effect has 
to be included, (c) The fine structure splittings may also 
be modulated by external electric and/or magnetic fields 
which are not discussed in this paper, (d) Our model 
can also be extended to investigate the initial optical spin 
polarizations J^ir-i^ 
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Appendix A: The Exciton Hamiltonian 

Here we write down the explicit form of the exciton 
Hamiltonian H C J},, in Eq. $2$ , which is recovered 

following the way laid out by Pikus and Bi r 41 ' 42 



1 2 ) — [Hm'mfel^n'n + H^, {k 2 )5 m : m 

ri r 2 J 

+ [7 ch (n - r 2 )6 m , m 6 n , n ] <5(n - r[)6(r 2 - r' 2 ) 
+ AIO„, T ' 2 ) , (Al) 



where k = — isj and 



tf oh (n - r 2 ) 



Atteqk n - r 2 



(A2) 



ff^ m (k!) = [E m (k ) + |^-k 2 



(kl • 7T m ' m ) 

m 



ft 2 (ki ■ 7r m ' n »)(ki ■ TT n " m ) 



in 



(A3) 



#n'n(k2 



E° m -E°„ 
H^ n@n ,(-k 2 ) = -H«, n (k 2 ). (A4) 



Here H^ n , m is in the usual form derived from the k • p 
method up to the second order. We have 7r = p + 
[tr x (vK))] with Vo standing for the lattice poten- 
tial. 7r m „ denotes the matrix element of 7r between the 
two Bloch functions V'mk ( r ) an d "0nk o ( r )- In this paper, 
lattice potential Vo is assumed to have spherical symme- 
try O stands for the time reversal operator. 

The e-h exchange interaction is decomposed into LR 
and SR parts: 



with 



1 1 *2 

ri r 2 



^LR ; / 1 1 1 2 

i'i r 2 



jjSR / I 1 1 1 2 

ri r 2 



(A5) 



ri r 2 



a/3 



On' 



d 2 

— ^—gUin-i^) 



x 5(n — r 2 )S(r[ 
r' r 1 . 



1 x 2 

ri r 2 

x 5(r 2 



~r' 2 ), (A6) 
<f(n - ra)5(n - r' x ) 

(A7) 
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In these equations, we have 



U(r) 



7T m'0n' 7r enm 



to 2 (2& -i5°,) ! 



47re K|r| ' 



(A8) 



(A9) 



and represents the eigenenergy of the s band at the 
point k = ko ; V is the volume of the bulk material which 
comes from the normalization of the Bloch function and 

x Otp n ' ko (r 1 )^ mko (r 2 )dr 1 dr2. (A10) 



Appendix B: The Bloch Functions 

The Bloch functions of the T§ conduction band take 
the formZ2 

\\\)c=\S)a, \\,~\) 

where a (/3) denotes the spin- up (down) state and \S) rep- 
resents the s-like conduction band Bloch function. The 
Bloch functions of the Tg and Ty valence bands are writ- 
ten asZ2 



(Bl) 



|^ = _-L(|x) + i|r))a 



2' 2 h 
3 1 
2'~2 
3 3 



1 

1 

i 



2' 2 7 



1 

71 



(B2) 

-(\X)+i\Y))/3 + 2\Z)a\, (B3) 
X)-i|Y»a + 2|Z)/3], (B4) 

(B5) 

X) + i\Y))(3 + \Z)a] , (B6) 



2 )„ = -j=QX) - i\Y))p, 



\-\)v = -^W *\ Y )) a - we] > ( B7 ) 

where |A), |Y) and |Z) are the p-like valence band Bloch 
functions which are real according to the phase conven- 
tion in accordance with the time-reversal symmetry. Af- 
ter taking the time reversal operation, we have 



3 3, 



1 



-,-)„ = - — QX)-i\Y))!3, 
l,l) v = l=[(\X)-i\Y))a + 2\Z)0\, 
\-\)v = ^=[(\X)+i\Y))P-2\Z)a], (BIO) 



(B8) 
(B9) 



l=(\X)+i\Y))a, 



3 3 

±,±) v = ±[(\X)-i\Y))a-\Z)0\, 
\-\)v = ~[(\X)+i\Y))0 + \Z)a]. 



(Bll) 
(B12) 
(B13) 



Appendix C: Construction of basis functions with 
direct Coulomb interaction explicitly included 

With the confinement, the diagonal part of the exci- 
ton Hamiltonian H c ^, n , f^J i.e., excluding the e-h 
exchange interaction, can be written into 



Hr> = Hp + Hh + Hi 



Coulomb 



H, 



confinement; 



(ci) 



where H e is the electron Hamiltonian in the form H e = 
h 2 k 2 /2m* e with to* being the effective mass of the con- 
duction electron and Hh is the hole Hamiltonian. From 
Eq. (fl4|) . we see that Hh is diagonal in the 4x4 matrix 
representation and the quasi-spins are not coupled. So in 
the following we omit the spin degree of freedom of both 
electron and hole. Hcouiomb and ^confinement are the di- 
rect Coulomb interaction and the confinement potential, 
given in Eq. (|A2|) and Eqs. ([3T |) -(|32 ]I . The eigen equation 
for the envelope basis function is constructed as 



H D \eh) = E\eh), 
(ri,r 2 |eh) = /(n,r 2 ). 



(C2) 
(C3) 



When a strong confinement is applied along the z- 
direction so that only the lowest electron/hole subband 
is relevant, one has 



/(n,r 2 ) = *(r lM ,r 2 ||)€(^i)C(^2), 



(C4) 



where £(z)[£(z)] = y^-sin(j £ ) stands for the electron 

(hole) envelope function in the z-direction. After multi- 
plying both sides of Eq. (|C2p with £(zi)£(z 2 ) and inte- 
grating over z\ and z 2 , one comes to 



(C5) 



with 



Hn = 



2ml 



2 

P 2 \\ 

2m* h 



i * / 2 2 . 2 2\ 



+ ?rn* h \\("lh4 +"l h yl) + V(r\\) + E% + £* (C6) 



in which 



to 



7i + 72 ' 

wo 
7i - 272 ' 



Et 



h 2 n 2 
2m*M 



2^2 



H 2 n 



2m* I 2 '' 



[X\ - £2,2/1 - 2/2), 



V{r{) = - / dzi 

l z JO JO 



dz 2 U ch (n - r 2 ) 
x {sm(irz 1 /l z )sin('KZ 2 /l z )} 2 . 



(C7) 
(C8) 

(C9) 



Here, to^m and m* h z are the effective masses of the heavy- 
hole in the plane and in the z-direction, respectively; 
E% and E h z are the subband energies resulting from the 
strong confinement in the z-direction. 
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After separating the coordinates of electron-hole pair 
in the plane into the center-of-mass and relative parts: 



r l| = I-! || - r 2 || 

m e r l|l+ m hj| r 2|| 



R 



K|| = ki|| +k 2 || 

_ ™h|| k l||-™e k 2|| , 



m t+ m U 



K. " rn t~\~ rn h\\ 

( C1 °) 

Ho is separated into two parts: Hp = Hq + H', where 



Hq 



h 2 l 
7T- ~{ k l + k y ) + -m jU w ave (a; + y 2 ) 

iLitib II £ 



1 „, 1 



1 



(K 2 x +K 2 ) + V( ril ), (Cll) 



2(m* +m*„) 



in which 



m 



m* 2 — m*m* h ,, + m 



CJnp /if 



TO e"*fc|| 



m* 2 — m*m* h n + mjj 



(C12) 
=**)~\ (C13) 
\ (C14) 



if' is further constituted of two parts: 
with 

b a 

m* m* n 

A [(W^ - Wave) 2 ' 2 + ( w yr - w avo)y 2 ] >(C17) 



(C15) 



= ™o(^r - ^x^x + <4>»n ( C16 ) 



Hl2 
where 



ml 



</ t || 



m* + m* 



m* + m A || 



J y0 



m l 2 ' 



(C18) 
(C19) 



Treating H' perturbatively, we see from Eq. (IC11I) that 
the center-of-mass and relative motions of the electron- 
hole pair are now decoupled and we are able to write the 
in-plane wave function as ^r^^r/jy) = ip(X, Y)cp(rn). 
Here tp is the eigenfunction of the 2D harmonic-oscillator 
potential 



(C20) 

where 



xH ny (p x Y), 

h 



h 



mtm* h J, 2 x 



-j Wye 



lm* e m* hli q 



(C21) 



(m* +m* H )uj xc 



, ft 



(m* e + m* H )u yc 
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FIG. 16: (Color online) The dependence of doublet splitting 
energy on the dot anisotropy in an GaAs QD with l z — 3.2 nm 
and Iq = 8.8 nm. r\ is set at 1 or to switch the direct 
Coulomb interaction on and off. The triangle points are taken 
from Fig. 2 in Ref. HI 



with its eigenvalue being E n ^ ny = (n x + \)Ttlo xc + (n y + 
Tj)hu) yc . A nx and A Uy are the normalization factors and 
H n (x) are the Hermit polynomials. 

According to Eq. (|C11|) . the relative part of the in- 
plane wave function can be expressed in polar coordinates 
as v>(ry) mn = -j=e lm ^ R n (p) with p = \r\\\. R n {p) is 
obtained by numerically solving the radial equation of 
the relative motion of the electron-hole pair in the real 
space. 

Overall, the basis function for Fm'n'( r i> r 2) m Eq. © 
is constructed as 



; (ri,r 2 ) = ii) nx n v (X,Y)ip mn (r^(zi)C{z 2 ), 

(C23) 

which is the eigenfunction of Hq with 



Hq — Hd — H . 



(C24) 



The exciton Hamiltonian H eh is diagonalized under the 
set of basis {\f mnn!D n y ) ® \ciVj)}, where \ciV ) acts as 
quasi-spin. In this way, the exciton eigenstates and eigen- 
values are obtained with the direct Coulomb and the ex- 
change interactions fully accounted. 



Appendix D: Comparison with results in Ref. 69 

We calculate the doublet splitting energies in single 
GaAs QD with rj = (the direct Coulomb interaction 
is switched off) and with rj — 1 (the direct Coulomb 
interaction is fully included). The length of dot major/ 
minor diameters are evaluated as l x 
In 



Iq\/1 + £, which in turn give u> x u 



h/y/T+l and 
wg,(l + f) and 
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LJyjj = + £) with = -^r-jj- Here □ stands 

for electron or hole, to^ denotes the corresponding in- 
plane effective mass and £ = (7 y — l x )/lx represents the 
dot anisotropy. This is consistent with that in Ref. [69|. 
The size parameters are chosen carefully to simulate the 
model employed in Ref. [6^ according to the characteristic 
energies induced by the confinement. In respect that the 
confinement potentials are chosen separately for electron 
and hole in Ref. [6^, we deduce from the the parameters 



therein two sets of size parameters. In the calculation, 
the size parameters are set between the corresponding 
two values. We choose l z — 3.2 nm and la = 8.8 nm. 
The dependence of the doublet splitting energy on the 
dot anisotropy for r/ = is shown in Fig. [16] (red curve 
with □) which is extremely close to that from Fig. 2 in 
Ref. [69| (green A). As we see, the absolute values of 
doublet splitting energy calculated with rj = 1 are much 
larger than those with r/ = 0. 
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